Scale-free energy dissipation and dynamic phase transition in stochastic sandpiles 
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We study numerically scaling properties of the distribution of cumulative energy dissipated in an 
avalanche and the dynamic phase transition in a stochastic directed cellular automaton [B. Tadic 
and D. Dhar, Phys. Rev. Lett. 79, 1519 (1997)] in d = 1 + 1 dimensions. In the critical steady state 
occurring for the probability of toppling p > p*=0. 70548, the dissipated energy distribution exhibits 
scaling behavior with new scaling exponents te and De for slope and cut-off energy, respectively, 
indicating that the sandpile surface is a fractal. In contrast to avalanche exponents, the energy 
exponents appear to be p-dependent in the region p* < p < 1, however the product {te — 1)De 
remains universal. We estimate the roughness exponent of the transverse section of the pile as 
X = 0.44 ± 0.04. Critical exponents characterizing the dynamic phase transition at p* are obtained 
by direct simulation and scaling analysis of the survival probability distribution and the average 
outflow current. The transition belongs to a new universality class with the critical exponents 

= 7 = 1.22±0.02, /3 = 0.56±0.02 and u± = 0.761±0.029, with apparent violation of hyperscaling. 
Generalized hyperscaling relation leads to P+13' — {d—l)iy±, where (3' = 0.195±0.012 is the exponent 
governed by the ultimate survival probability. 



I. INTRODUCTION 

Open driven systems exhibiting self-organized critical 
states are usually modeled by sandpile-type automata, 
in which sand grains are added slowly to the system and 
its evolution is monitored in terms of collective sandslides 
(avalanches). With additional complexity due to stochas- 
tic character of the relaxation rules , stochastic cellu- 
lar automata models have proven useful in understanding 
certain aspects of granular flow |^ |^ in realistic granular 
materials and in stochastic biological processes [|. On 
the other hand, sandpile automata models are interesting 
from the theoretical point of view since both role of the 
dynamic conservation law — conservation of the number 
of grains — and emergent spatial structures can be easily 
monitored. 

Recently a model with probabilistic toppling and di- 
rected mass flow has been proposed [Q, in which the 
probability of toppling p represents a control parame- 
ter originating either from random variations of stick- 
ing properties between grains, or from stochastic exter- 
nal conditions (wetting and drying properties). Another 
realization is related to stochastic processes in biological 
systems such as random dispersion of particles, which are 
added from the outside and evacuated from the system 
only when its response lasts longer than a fixed time Tq, 
measured on the internal time scale of the process. Parti- 
cles are transferred among communicating cells according 
to probabilistic rules, however, for response times shorter 
than To they are held inside the system. Therefore, each 
cell contains a certain number of particles, which varies 
with time. The probabilistic character of the particle 
transfer between connected cells can be attributed to 



mechanisms which depend on general condition of the 
system. We consider the case Tq = L, where L is the 
linear system size. 

The relaxation rules of the model in d = 1 -f 1 dimen- 
sions are: if h(i,j) > he = 2 then with the probability p 
the relaxation occurs as follows: 

h{t,j) ^ h{i,j) - 2; h{i + 1, j±) ^ h{i + 1, j±) + 1. (1) 

Here h{i,j) is the dynamic variable, i.e., the height (num- 
ber of particles) at site (j,i), and (i -t- 1, j±) are neigh- 
boring downstream sites on a 2-dimensional square lat- 
tice oriented downwards. Due to directed mass transport 
the dynamics of this model is anisotropic leading to self- 
affine relaxation clusters in d = 2 (see Fig. 1). Thus the 
model can also be viewed as directional lines in d = 1-1-1 
dimensions, in which the instability can propagate in one 
spatial dimension back and forth, whereas the temporal 
dimension is strictly directed. The system is driven by 
adding particles from the outside, one at a time at a ran- 
dom position along the first (top) row. Only sites which 
are connected to toppled sites at the previous time step 
are considered as candidates for toppling. Perturbation 
can then be transferred from an active (toppled) site to 
two forward neighbors. Periodic transverse boundaries 
are assumed and all candidate sites are updated in par- 
allel. The probabilistic character of the relaxation rules 
produces a ragged structure of heights. [A transverse sec- 
tion through the pile is shown in Fig. 1 (top)]. Two top 
particles from the surface h{i,j) are taken away from the 
system when the active site is on the lower boundary, i.e., 
i = L. 

It should be stressed that according to the above re- 
laxation rules transport of grains is independent of the 
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relative heights of neighboring sites. Therefore, it may 
occur that at some sites height difference in the direc- 
tion of transport is negative, and thus the system per- 
forms work in order to maintain the transport. A nice 
example of this dynamic rules was found recently in bi- 
ological transport processes which are mediated by so 
called molecular motors Q . These are protein molecules 
that can use excess energy from the chemical reactions in 
the fuel (adenosine triphosphate) and perform mechani- 
cal work. 

It has been understood that the automaton exhibits 
self-organized criticality for the range of values of the con- 
trol parameter p > p*, where p* = p^^^ = 0.7054853(5) 
is the percolation threshold for the site directed percola- 
tion on the square lattice. Due to the dynamic con- 
servation law — conserved number of grains in the interior 
of the pile — and the probabilistic relaxation rules which 
are locally like the generalized site-bond directed perco- 
lation (see below), the system exhibits emergent spatial 
structure, as discussed in detail in Ref. [Q. Moreover, 
the avalanche exponents for the integrated probability 
distributions of duration P{T) ~ y-(n-i) g^j^j gj^g of 
avalanches, D{s) ~ s~^'^''~^\ are expressible in terms of 
the standard directed percolation exponents in d = 1 + 1 
dimensions as follows [0: 

Tt 1 = a - (d - IKdp - (/J/^^iOdp ' (2) 

and Tg = 2 — 1/rt, and the anisotropy exponent ^ (^up. 
Subscript DP refers to directed percolation, and /3 and 
I'll are the critical exponents for the order parameter and 
parallel correlation length, respectively. 

In the present work we extend the study of the model of 
Ref. 01 in two ways: (1) We study the probability distri- 
bution of the potential energy dissipated in a relaxation 
event (avalanche) due to grains drop from higher to lower 
positions at the fluctuating sandpile surface. This distri- 
bution is unique for the dynamic sandpile models with 
ragged spatial structure and has no counterpart in the 
directed percolation processes. Thus we expect that the 
exponents characterizing its scaling properties are new. 
(2) We analyze the behavior of the system close to the 
phase transition point in terms of (i) scaling properties 
of the survival probability distribution for p < p* and 
(ii) by determining the time averaged outflow current 

< J{p) >, which behaves as an order parameter. The 
outflow current results from avalanches which last longer 
than the system size T > L. The average is taken over ex- 
ternal time scale, which is measured in number of added 
grains. The internal current at time T < i is defined as 

< j{T,p) >~ T^°'m{T,p), where ■m{T,p) is the average 
flux of particles at time T, and is the probability 
that an avalanche survives T steps. In the steady state 
outflow current < J{p) > balances the input current, 
which is one particle per time step, and thus it is equal 
to one. For p < p* the system ceases to conduct parti- 



cles, and the outflow current drops to zero for £ ^ ^{p) 
and finite lattice size L as 

< J{p) >~ {Sprg{L/ap)) , (3) 

where Sp = (p — p*)/p* measures the distance from the 
steady state and ^{p) ~ ((5p)~''ii is the parallel correla- 
tion length. The set of critical exponents is determined 
by the appropriate scaling fits and using the scaling rela- 
tions that are valid in the present dynamical model (see 
Sections IV and V). 

The organization of the paper is as follows: In Sec. 
II the phase diagram of the system is calculated nu- 
merically for finite lattice size L. In Sec. Ill we study 
scaling properties of the distribution of energy dissipated 
in avalanches. In Sec. IV we present detailed numerical 
analysis of the phase transition. Summary of the univer- 
sal scaling exponents and the discussion of the results is 
given in Sec. V. 

II. PHASE DIAGRAM 

Due to probabilistic dynamic rules the avalanches in 
this model show a fractal structure; an example is shown 
in Fig. 1 (bottom). In the limit p = 1 this model re- 
duces to the deterministic directed model with compact 
avalanches, which has been introduced and solved exactly 
in Ref. 

As discussed in detail in Ref. , for < p < 1 the re- 
laxation rules at each site of the system may be visualized 
as the rules of a Domany-Kinzel cellular automaton |l^] 
of generalized site-bond directed percolation (DP), with 
probabilities Pi and P2 that a toppling occurs if one or 
two particles, respectively, drop at that site. According 
to Eq. (0), we have Pi = pp, where p is the probability 
that the site has height h > 1, and P2 = p hy definition. 
In contrast to DP, in the present dynamic model the state 
is being systematically built up after each avalanche, and 
the probability p was found to vary with the distance £ 
from the top of pile as p{£,p) = p*{p) — A{p)£~^ . 
Here x = ^/v^^ (see inset to Fig. 2) is the inverse par- 
allel correlation length exponent of the directed percola- 
tion jl^. In the above formula p*{p) are the values of 
p{£,p) reached ai £ ^ 00. Notice that the distance £ 
from the top of the pile is equivalent to the duration T of 
avalanches. In Fig. 2 we show time-averaged < p*{p) > 
vs. p obtained numerically and averaged over lower third 
of the lattice with L = 100, for various values of p. Two 
types of initial conditions are used: (a) full lattice (all 
lattice sites are occupied by at least one particle), and 
(6) half-full lattice (half of the sites which are selected 
randomly have zero heights and the rest of sites are oc- 
cupied). In both cases, for the initial set of probabilities 
(p, p*) in the region to the right of p* the system self- 
organizes (after some transient time) to sitting close to 
the DP critical line (cf. Fig. 2). Left of the line p = p* 
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different initial conditions lead to separate final states. 
We used 8 x 10^ time steps for each point. Notice that 
due to finite size of the lattice < p*{p*) > is still some- 
what smaller than one (indicated by dotted line in Fig. 
2), and that a dynamical hysteresis occurs in the region 
0.5 < p < p* |l3|. At p* an instability — building up of 
heights — starts at lower boundary of the pile and pro- 
liferates inside, reaching the first row for probability of 
toppling exactly p = 1/2. Apart from the finite-size ef- 
fects, the phase diagram in Fig. 2 is in agreement with 
general theoretical considerations given in Ref. [Q. 

III. DISSIPATED ENERGY DISTRIBUTION 

Due to the probabilistic character of the dynamic re- 
laxation rules in Eq. (|l|) for p < 1 and the conserva- 
tion of particles in the interior of the pile, our dynamic 
model exhibits the emergent spatial structure |Q, which 
is characterized by rough surface h{i,j) embeded in 3- 
dimensional space. Therefore, mass transport in the pre- 
ferred direction takes place along a rough surface. Con- 
ditions for the potential energy dissipation are fulfilled 
when the height difference 9||/i± = h{i, j) — h{i + l, j±) — l 
along the direction of transport is positive. More pre- 
cisely, since relaxation at a site involves two particles, the 
energy is dissipated at that site when the sum E(z,j) = 
(9||/i_ - 1 > 0. 

Two comments are in order at this point: (1) Although 
the driving force of the grain transport in this model is 
not a gradient of the potential energy, as discussed in the 
Introduction, we believe that cumulative potential en- 
ergy dissipated in an avalanche is an interesting quantity 
which is directly related to temporal fluctuations of the 
emergent structure in the real space. We concentrate on 
the properties of the energy distribution near the phase 
transition point p = p*, where the surface exhibits dra- 
matic fluctuations, and only briefly discuss scaling behav- 
ior for p > p*. (2) We distinguish between energy dissi- 
pated at a fixed site, through which various avalanches 
run, and the energy dissipated at different points in the 
whole avalanche. Here we study scaling properties of the 
latter quantity. Notice that the condition E{i,j) > is 
fulfilled at a set of points S which is a random subset 
(see Fig. 1) of the avalanche size s jwj. Total energy 
dissipated in an avalanche is then E = E{i,j). 

The probability distribution of the dissipated energy 
P{E) is found to obey a power-law behavior with the ex- 
ponent te for p* < p < 1 and the following scaling form 
is satisfled 

P{E, L) = L^'-'^-^^^'^ViEL-^") (4) 

with (te — 1)-D_e = Tt~^, where Tt — 1 = a is the survival 
probability distribution exponent. In Fig. 3 we show the 
integrated distribution of dissipated energies for p — p* 



and for four different values of lattice size L. The slope 
gives the exponent te ~1 = 0.24, and the finite-size scal- 
ing plot according to (|4|), which is shown in the inset 
to Fig. 3, is obtained with a =0.45 and De =1.84. In 
contrast to the survival probability distribution and size 
distribution of avalanches, the distribution of dissipated 
energy cannot be defined in the directed percolation pro- 
cesses, and thus the exponents te and De are new and 
are not directly related to the DP exponents. Moreover, 
we find that the exponents te and De are p-dependent, 
for instance, for p =0.8 we obtain te = 1.27 and De = 
1.66, and te = 1-29 and De = 1.55 for p =0.9. However, 
a combination of these exponents can be related to the 
survival probability exponent a via the scaling relation 
{te — 1)De = a, which holds in the SOC states, where 
a is the universal exponent expressible in terms of DP 
exponents via Eq. (H). This scaling relation is satisfied 
within numerical error bars (estimated as ±0.02) for all 
values of p in the region p* < p < 1- In the limit p =1 
the critical state is exactly known ||l^ and consists only 
of the heights h = 1 and h = Q. Consequently, dissipated 
energy is bounded to values flfe integer values, and the 
distribution P{E) has the same scaling exponents as the 
size of avalanches distribution. 

For p* < p < 1 emergent spatial structure appears 
due to both stochastic dynamics and the conservation of 
number of particles in the interior of the pile ^ . At the 
edge of the scaling region {p = p*) we find that the aver- 
age height (averaged in the transverse direction) increases 
with the distance £ from the top row as < h{£) >w a£^, 
with b = 0.59±0.02 (cf. Fig. 1 (top)). Therefore, for large 
£ there is a finite probability of large heights, however, 
at the same distance £ some site have height zero, since 
the system is in the stationary critical state (in the op- 
posite the avalanche would propagate as a directed per- 
colation cluster, which violates stationarity condition). 
Thus, 9|| h also increases with £ and becomes unbounded 
for £ oo. Thus, the energy cutoff has additional non- 
trivial ^-dependence, which is not contained in the £- 
dependence of the avalanche size cutoff, indicating that 
the sandpile surface is a fractal at p — p*. In the in- 
terior of the scaling region {p* < p < 1), the average 
height remains finite and not a function of £, however, 
height distribution does strongly depend on p. We flnd 
that the width of the height distribution w{p) increases 
smoothly with decreasing p from w = latp=ltoa 
fast diverging function a,t p —>■ p* . Neighboring sites are 
weakly correlated since the dynamics is governed by the 
critical height rule only, and thus dissipated energy at 
a site E{i,j) also depends on p. We checked by direct 
calculation that the distribution of dissipated energy at a 
flxed site in the interior of the pile, P{Esite), taken over 
2 X 10^ avalanches exhibits strong p-dependence. It is 
an exponential function of width Wes , which is increasing 
smoothly with decreasing p and becomes nearly power- 
law a.t p = p*. We believe that the p-dependence of the 
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height distribution is the origin of the observed nonuni- 
versaUty of the energy exponents. On the other hand, size 
and duration of avalanches are governed by the probabil- 
ities p and p*, which sit always at the DK critical line 
in Fig. 2. It should be noticed that the probability p{£) 
does not depend on particular values of heights h > 1, 
and thus the avalanche exponents remain universal (cf. 
Eq. (D). 

We find that the distribution of mechanical work done 
by the system exhibits a curvature and not a power-law 
behavior. 

IV. DYNAMIC PHASE TRANSITION 

In the region below p* the system ceases to conduct 
and starts accumulating particles. As a consequence the 
critical steady state is lost (see detailed discussion in Ref. 

) and the probability distributions exhibit exponential 
cut-offs with finite correlation length, depending on the 
distance from p*. In Fig. 4 we show the survival prob- 
ability distribution P{T,p) for few values oi p < p* and 
L =200. In general, a distribution P{X,p,L) satisfies 
the following scaling form in the subcritical region 

P{X,p,L) = , (5) 

where 6p = {p* ^p) /p* and X stand for T, s, or i?, respec- 
tively, and Dx is the corresponding fractal dimension. In 
the case of distribution of durations P{T,p,L) we have 
Dt = z the dynamic exponent, and z = 1 in the present 
model. Therefore there are no finite-size effects in the 
survival probability distribution, which makes it partic- 
ularly suitable for the subcritical scaling analysis. In the 
case of size and energy distributions, one is restricted to 
values of p and L such that the condition (Jp)""^!! /L 
is satisfied. In the inset to Fig. 4 the scaling collapse ac- 
cording to Eq. (^) of the survival probability distribution 
is shown, where we have used a = 0.45 and zv\\ — 1.22. 

Another way to study the dynamic phase transition is 
by direct measurements of the order parameter, i.e., the 
time-averaged outflow current < J{p) >. In the critical 
steady state < J{p) >= 1, thus balancing the average 
input current. Below the transition point this balance is 
lost. The outflow current decreases reaching zero at some 
lower value of p, which depends on the system size L [see 
Fig. 5 (bottom)]. For different lattice sizes we expect the 
following scaling form to hold: 

< J{p, L) >= L-fZ-^ii J(Li/'^,i _ p*yp*) . (6) 

This scaling form follows from general scaling proper- 
ties of the internal current for T < L, that reads: 
< j(T,p,L) >~ L^^jiL^^-^iiip - p*)/p*,L-^T). By 
choosing L ^ [p/p* — 1) ^^i' = ^ and having defined 
the exponent /3 in Eq. (^, we find that the anoma- 
lous dimension \j = —l3/v\\. Therefore < j(T,p,L) >— 



L-'^/^ii g(Li/^ii {p-p*) /p*, L-^T). In the stationary state 
for p > p* , correlation length ^ — > cx) and the first 
argument in Q can be neglected. For T we ex- 

pect that the scaling function Q behaves as a power, i.e., 
< j{T,p,L) >~ const X T-"L^"-'5/''ii , which should be 
independent on L, thus leading to = za. For the 

outflow current, however, we have T > L and second ar- 
gument of Q can be neglected. Then for p < p* one gets 
the expression (|). Taking L^'^ii _ {p-p*)/p* leads to 
Eq. 1). 

The scaling plot according to Eq. (|6|) is shown in Fig. 
5 (top) where we have = 0.45 and = 0.83. 

Together with the above results and observing the er- 
ror bars for a from we estimate the exponents as 

= 1.22 ± 0.02 and /3 = 0.56 ± 0.02. 

Notice that the above values of the exponents z/|| and /? 
are close to the values 1.28±0.06 and 0.58±0.06, respec- 
tively, obtained by Monte Carlo simulations for d = 3 
dimensional directed percolation in Ref. |l5|. The rea- 
son for this similarity lies in the altered character of the 
dynamics of our model below the transition. Namely, 
the probability p* of having height h > 1 reaches unity 
(in the limit of large L) at p* . The consequences of this 
are twofold: (i) the threshold character of the dynam- 
ics is lost ioi p < p* , i.e., each site in the lattice satisfies 
the threshold condition h > he = 2 when a single particle 
drops on that site; {ii) since p — 1 we have PI = P2 = p, 
thus the system spreads the perturbation in the (i?|| , R±) 
plane with the probability p and builds up heights with 
probability q = l—p- Therefore, for 0.5 < p < p* we have 
a dynamic model in which an avalanche propagates effec- 
tively as a cluster in a 3-dimensional directed percolation 
with finite widths ^|| , in the plane, and percolating in 
the vertical direction. However, there are considerable 
differences between these processes and conventional 3- 
dimensional DP, leading to generally different set of ex- 
ponents, as discussed below. An added particle moves 
along rough surface h{R\\,R±), which fluctuates inside 
the correlated region (^||,^_l). The internal time scale 
becomes bounded by finite ^||, and the system percolates 
for t —f oo, where t is now the external time scale (mea- 
sured by the number of added particles). Flights of par- 
ticles along the rough surface are proportional to local 
height gradients which are usually larger than one, 
in contrast to contact percolation processes. However, 
the average height of the pile < h > grows exactly by 
one unit with each added particle as a consequence of 
the conservation of number of particles. 



V. DISCUSSION AND CONCLUSIONS 

The dynamic model with stochastic relaxation rules 
of Eq. (|^) exhibits the universal self-organized critical- 
ity for the range of toppling probabilities p < p* < 1, 
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and the dynamic phase transition at p = p*. As dis- 
cussed in detail in Ref. the avalanche exponents a 
for survival probability, r = r, — 1 for integrated clus- 
ter size distribution, and ^ for the average transverse 
extent of clusters, are expressible in term of standard 
directed percolation exponents in all dimensions via Eq. 

. Here we have shown that the dissipated energy dis- 
tribution, which is peculiar to the dynamic model and 
has no analogue in the directed percolation processes, 
is described by a new exponent te and corresponding 
fractal dimension De- These exponents appear to be p- 
dependent, however, their product is related to the uni- 
versal survival probability exponent due to scaling rela- 
tion (te — 1)De — OL. Power-law behavior of the distribu- 
tion D{E) a,t p = p* indicates that the sandpile surface 
h{i,j) is a fractal. We estimate the roughness exponent x 
by measuring the time averaged height vs. transverse di- 
mension of the pile < h{Rtr) >^ Rfr at various distances 
£ from the first row [see Fig. 1 (top)]. By box counting 
we find that the contour curve of the perpendicular sec- 
tion through the pile for p ^ p* has fractal dimension 
df = 1.44 ± 0.045, leading to x = df -1 = 0.44 ± 0.04. 
Error bars are estimated from several separate measure- 
ments at different sections. The roughness exponent ap- 
pears to be larger than the one measured in ricepile model 
with critical slope rules, where it was found xrp =0.23 
1^,^. In the steady state the height fluctuates around 
the average value < h >= 19±6, increasing by one 
unit at boundary sites of an avalanche and decreasing 
by one unit at transport sites with one preceding active 
neighbor. The flame-like profile in Fig. 1 (top) indicates 
individual site fluctuations, in agreement with the criti- 
cal height rules in our model. Closeness of the exponents 
X « a indicates that the number of sites at which the pile 
grows is on the average equal to the number of transport 
sites, i.e., the avalanches have almost no compact parts 
(cf. Fig. 1). Below the transition point the pile grows 
indefinitly for f ^ oo, as discussed in Sec. IV. 

The dynamic phase transition at p* is characterized by 
the exponents of the order parameter /3, and the paral- 
lel correlation length z/|| , the numerical values of which 
appear to be very close to those of the directed perco- 
lation in d -|- 1 dimensions. However, the exponent 7 
for the order-parameter fluctuations and the exponent k, 
which describes the average cluster growth at time T as 
m{T) ~ T'^, appear to be different from d = 3 DP ex- 
ponents (see Ref. [|5[). A complete set of exponents is 
given in Table I. 

With regard to the exponents in Table I we would like 
to point out the following: (1) In the critical steady state 
the average cluster growth balances the average input of 
particles, i.e., ■m{T) ~ 1, leading to k = 0. The scaling 
relation k = 7/1^11 — 1 = D|| — 1 — a holds, thus we have 
7 = ; (2) At the dynamic phase transition the hyper- 
scaling relation (HS) 2(3 + j = [1 + {d — 1)C]'^\\ appears 
to be violated , in contrast to the thermodynamic 



DP phase transition (cf. Table I). The avalanche expo- 
nents are determined in Ref. 0. The exponents for DP 
in d = 2 are taken from Ref. |9| and the correspond- 
ing avalanche exponents tjjp and Djjp are calculated 
using the scaling relations. (Notice that due to the pres- 
ence of anisotropy, at least three exponents should be 
known in order to determine completely the universal- 
ity class.) The following scaling relations are valid both 
in DP and in dynamic models: — Tt ^ I = a; 

P/{f3 + 1) = Ts - 1 = T- C = v±/v\\, and D\\v\\ = /? -h 7. 
The hyperscaling violation exponent is defined via 
/3 -I- 7 = [1 + (d - 1)C - + , which together with 
the above equations leads to fl = D^^ — 1 + a — {d — 1)^. 
Here D^^ is the fractal dimension of the size of relaxation 
clusters measured with respect to the length parallel to 
the transport direction. Using the fact that D[| = 1 + a 
in the dynamic model, we may write f2 = 2a — (d — 1)^. 

In the directed dynamic processes it is useful to define 
another exponent /?', such that the generalized hyper- 
scaling relation |p^-|20|| 

a{l + f3'/f3) + {d- 1)C , (7) 

is satisfied. Here /?' is related to the ultimate survival 
probability (the survival probability of a cluster grown 
from a fixed seed), whereas P governs the usual order- 
parameter — stationary density of active sites. Recently 
exponential inequality (3' ^ (3 was found in models with 
multiple absorbing configurations jl^ and in branching 
annihilating random walks with even parity | |20| |. Our 
results in this paper suggest that in the self-organized 
dynamic critical states in directed models k = and 
/3 7^ /?' is always satisfied. We have 

^ A _ > . (8) 

Using the above scaling relations and Eq. (||) we find that 
(3'/v\\ = aDP=0.159 (cf. Table I). Therefore, the HS vi- 
olation exponent f2 is given by the difference between 
the survival probability exponents in the dynamic model 
and the underlying directed percolation as 17 = a — cuDp^ 
which turns to be equal to the cluster growth exponent 
of DP, i.e., = KDp. One can also define a new con- 
figuration exponent 7' via ^' = 1 + such that 
the relation 2/3' + 7' = (1 + (d — l)C)i'|| is satisfied. By 
inserting the above expression for f2 into 7', we have 
7'/i^j| = 1-1- Kp,p = 1.314. Using value of from Table 
1 wc find 7' = 1.586 ± 0.026, and (3' = 0.195 ± 0.012. 
The origin of the difference between a and aup lies in 
the dynamic conservation law, which leads to the emer- 
gent spatial structure and dependence of the branching 
probability p{tj on distance as discussed in Ref. [Q. 
Our present results suggest that the dynamic conserva- 
tion law is also responsible for the new universality class 
of the dynamic phase transition at the edge of the critical 
region. 
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TABLE I. Critical exponents in stochastic directed model 
(PH) and in directed percolation (DP) in d = 2 dimensions. 
Also shown are exact exponents in p = 1 limit (DR). 
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DP 


0.159 
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1.472 


0.634 


0.314 


0.276 


1.734 


2.278 
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FIG. 1. Bottom: Avalanche in 1 4- 1 stochastic model 
for p = p* and L =128. Transport direction is from left to 
right. Dark points represent sites at which system performs 
work, and at light dark points it dissipates energy. Top: (a) 
Transverse section of the pile at distance ^ = 96 from the 
top. Heights are averaged over 5 x 10*' runs, (b) Height (av- 
eraged in transferse direction) vs. parallel distance from the 
top in the stationary state. Dotted line: 1.6^°'^®. Heights are 
measured in number of grains and R is given in lattice sites. 
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FIG. 2. Phase diagram in the space of probabilities (p, p*) 
obtained numerically for a lattice with L = 100. Open sym- 
bols correspond to initial conditions along the line (p, 1), and 
filled symbols to initial conditions along the line (p, 1/2). Dot- 
ted line indicates value < > at p*. Inset: p{(.) vs. £ for 
(top to bottom) p — p* , 0.8, and 0.9, with fit lines described 
in the text. 
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FIG. 3. Double logarithmic plot of the integrated distri- 
bution of dissipated energies D{E, L) vs. E for p = p* and 
for L = 12,24,48, and 96 (left to right). Inset: Scaling plot 
according to Eq. (0) with a = 0.45 and De = 1-84. 
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FIG. 4. Double logarithmilP jJlot of the integrated distribu- 
tion of durations P(T,p) vs. T for L=200, and p=0.695, 0.69, 
0.68, 0.67, and 0.66 (top to bottom) in the subcritical region. 
Inset; Scaling plot according to Eq. (^ with Sp = (p* ~ p)/p* 
and exponents a = 0.45 and zvu = 1.22. 
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FIG. 5. Average outfiow current < J(p) > vs. p for vari- 
ous lattice sizes (bottom). Finite size scaling plot according 
to Eq. (|) with /3/i/|| = 0.45 and l/i/|| = 0.83 (top). 
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